Formulae for partial widths derived from the Lindblad equation 
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A method for calculating partial widths of auto-ionizing states is proposed. It combines either a complex 
absorbing potential or exterior complex scaling with the Lindblad equation. The corresponding classical rate 
equations are reproduced, and the trace conservation inherent in the Lindblad equation ensures that the partial 
widths sum up to the total width of the initial auto-ionizing state. 
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I. INTRODUCTION 

Resonances, or meta-stable states, play a crucial role in several physical phenomena - ranging from solid state physics ||1 



Resonances, or meta-stable states, play a crucial role in several physical phenomena - ranging from solid state physics 1 1-4 J 
£^ and plasma physics Jl,01 via atomic and molecular physics 0-03] to nuclear physics |[Illll2tl . The population of such states 
are frequently assumed to follow an exponential decay law. However, we know that this "law" is broken - both at short and 
long time scales. The former is due to the quantum Zeno paradox, whereas the latter is related to the fact that there exists a 
£f*^ ' lower threshold for the energy of the system lfl3tl . In a large intermediate time region, however, exponential decay remains a 
£SJ , perfectly valid approximation. In such a context it is customary to assume that as the population of some initial meta-stable state 
decreases exponentially, also the population of the (possibly) various decay products follow rate equations. When making this 
assumption certain interference effects are neglected, and one may say that this is a semi-classical assumption. This suggests 
that an approach based on a master equation describing a density matrix rather than a pure state is appropriate lfl4l - [T6ll . 

Resonance states may appear as solutions of the time-independent Schrodinger equation if we lift the restriction that the 
solutions should be square integrable and impose outgoing boundary conditions. Asymptotically, these outgoing waves are 
characterized by complex wave numbers. Consequently, also the eigen energy becomes complex. This is not in violation 
Z-> of the Hermicity of the Hamiltonian since Hermicity only applies to the space consisting of normalizable wave functions; 
O" 1 these complex energy solutions diverge. This is obviously a rather undesirable feature if we want to consider resonance states 
\ explicitly. To this end, methods which involve explicitly non-Hermitian terms in the Hamiltonian have proven very useful. When 
■ it comes to resonances, the usefulness lies in the fact that such terms are able to "kill off" the divergence such that also resonance 
£> ' states may be represented by normalizable functions. Such explicit non-Hermicity may be introduced in various ways. It may, 
e.g., arise by introducing an artificial imaginary potential which is zero in some interior region and increasing towards some 
boundary (a complex absorbing potential - CAP) Il7l - [l9tl . Other frequently used techniques are (uniform) complex scaling 
J20I - I22I1 . exterior complex scaling (ECS) l23ll24ll and smooth exterior complex scaling ll25il26ll . 

In literature, several examples of combining formalism for open quantum systems with explicitly non-Hermitian Hamiltonians 
are found, see, e.g., J27l - [32ll . The term "open" in regard to a quantum system indicate that the system is interacting with some 
environment. As the environment typically has a large number of degrees of freedom, a full solution of the total composite system 
can usually not be obtained. In such a context, it may be desirable to express the evolution of the smaller system alone where the 
effect of the environment on the smaller system is included somehow. This is typically achieved by invoking the Born-Markov 
approximation lO . The resulting equation of motion, which provides the dynamics of the reduced density matrix of the system, 
is frequently used in order to describe, e .g., relaxation due to the interaction with a radiation field J^,[13l, a damped harmonic 
oscillator l35tl or quantum transport l36l l37| |. Moreover, applications in which the system's own continuum is considered to be 
an environment are also found in literature I27ll3lll38tl . (The terms "open" and "unbound" are sometimes used interchangeably.) 
In lT28l l29ll it was demonstrated how the Lindblad equation may be used to describe the dynamics of spontaneously decaying 
particles, i.e. to systems in which the number of particles is not constant. It has proven useful to combine these ideas with the 
formalism of second quantization l30ll36ll . Master equations have also proven useful to incorporate the process of spontaneous 
decay in plasma physics; by introducing the relevant rates into the master equation, non-unitary processes involving auto-ionizing 
states may treated on equal footing with unitary processes involving radiation and collisions |5|,|g]. 

This work aims to combine non-Hermitian quantum mechanics with a master equation in order to describe the population 
dynamics of a system undergoing spontaneous decay to a system with fewer particles. We will make use of the ideas presented 
in ll30ll . which demonstrates how particle loss due to a CAP may be described in a consistent manner. This formalism allows us 
to disregard escaping particles and focus upon whatever is left. (In general this cannot be done using the Schrodinger equation.) 
Thus, we should be able to maintain information about the part of the system that remains bound. It is argued that the Lindblad 
equation is the proper framework for achieving this. Thus, we take the Lindblad equation to be our starting point; the equation 
of motion is not derived from any microscopical arguments. The remainder of the system after absorption is in general described 
by a density matrix rather than a wave function; as the degrees of freedom corresponding to escaping particles are "integrated 
out", some coherence effects are lost, see also e.g. |[39tl . 

We will assume that the system starts out in a resonance state and decay by emitting a single particle. Of course it may 
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rightfully be asked what is meant by "starting out in a resonance state" as such states are not physical. We will not enter into this 
discussion here, however, but rather simply exploit the fact that the resonances may be represented as square integrable eigen 
states of some (artificial) non-Hermitian Hamiltonian. We will assume that after emission of a particle, several final states are 
accessible. In this context, only Hamiltonians with no time dependence are considered. As it turns out, all time dependence is 
seen to follow the exponential decay law, and the resulting formulae for partial widths are time independent. 

The paper is organized as follows: In Sec. HI] the formalism for a CAP is presented. For completeness, the results of lf30ll are 
briefly outlined. In Sec. Ill Al it is demonstrated how the Lindblad equation reproduces the corresponding classical rate equations, 
and from this correspondence the partial widths are identified. The loss of coherence is briefly discussed in Sec. lIIBI In Sec. Hill 
the formalism is generalized so that it may also be applied to obtain partial widths using exterior complex scaling (ECS). In 
Sec.[IV]the applicability of the proposed formulae is discussed, and some remaining challenges are pointed out. Conclusions are 
drawn in Sec. [V] 



II. PARTICLE LOSS DUE TO A COMPLEX ABSORBING POTENTIAL 



Any non-Hermitian Hamiltonian may be written as the sum of a Hermitian and an anti-Hermitian part; 

H = H h -iH ah with H h = (H h y and H ah = (if ah ) t . (1) 

We will assume that all eigenvalues of H have non-positive imaginary parts. A sufficient condition for this is that iJ ah be positive 
semi-definite ll40l l4lll . In l30ll the anti-Hermitian part consisted in a CAP. It was demonstrated, expressing all interactions in 
terms of second quantization, that the Lindblad equation, contrary to the Schrodinger equation, is able to restore the dynamics 
of the remaining particles after absorption of other particles in a consistent manner. 

Suppose an A^-particle system is described on a numerical grid {.t;}, where X{ refers to the position, and also other degrees 
of freedom such as the spin of a particle. Here, a countable representation has been chosen for notational simplicity. With 
creation and annihilation operators cj and Ci, where cj creates a particle at "position" Xi and Cj annihilates a particle at x,, the 
Hamiltonian may be expressed in a manner which does not depend on the number of particles at hand: 

H = hkjclci + i ^ Vp q ^ rs clclc s C r . (2) 

kl pqrs 

We have here assumed at most two-particle interactions. The field operators obey the usual anti-commutation relations 

= 0, {ct,ct} = 0, {c i ,c]} = S iJ . (3) 
The CAP, which is a local imaginary potential — iT(x) with T > 0, may be represented analogously by 

f = ]Tr(x l )cJc l . (4) 

i 

Since the process we wish to describe should preserve total probability (trace), Trp = 1 Vi, and positivity, p > Vt, and the 
absorption process should be Markovian, our equation of motion should be on Lindblad form ll42j|43ll : 

ihp = [H, p] - i 7*,J {{A\Ai , P} - 2A lP A\) . (5) 

kl 

By comparison with the von Neumann equation for a non-Hermitian Hamiltonain, 

ihp = Hp- P H^ = {H h ,p]-i{H ah ,p} , (6) 

it is seen that H in Eq. (O should be identified with the Hermitian part of the Hamiltonian, H h , and the Lindblad operators 
should fulfill 

^ 7M 4A i= ff ah =f . (7) 

kl 

This leads to the equation 

iHp=[H,p]-i{f,p} + 2iJ2 r (%k)c k pcl (8) 

k 
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for the density matrix of the entire system. In this Fock space description, the total density matrix p does not correspond to 
a fixed number of particles. For an initial state with a well-defined particle number N, the total density matrix remains block 
diagonal, where each block corresponds to a sub-system consisting of n < N particles. The n-particle sub-system evolves 
according to 

ihp n = [H, p n ] - i jf , pnj + ihS[p n+1 ] , (9) 

where the source term 

2 x - t 

S[Pn+l\ = T / ,T(x k ) Cfc p n+1 c\ . (10) 
fc 

For the initial A^-particle sub-system, the von Neumann equation is reproduced with the "effective Hamiltonian" Hqap = 
H - it, 

ihp N = H CAP p N - p N i?cAP ' (11) 

which is equivalent to the Schrodinger equation if the initial state is a pure state, p(t = 0) = vpW^iJ/W |. 

Of course, the field operators c[ ' could also refer to single particle orbitals, Xi( x )' rather than points on a grid. In such a 
representation, the CAP is in general not diagonal, and Eq. (TTOt is rewritten as 

S[p n+1 ] = -^2(xk\F\xi) ci p n+ i c{ . (12) 

kl 

We would like to stress that this formalism is based solely on the Lindblad equation. The fact that we wish to describe 
a Markovian process in a manner which preserves trace and positivity, cf. j42l 14311 . is enough to justify the generic form of 
Eq. (0. We do not make any explicit reference to any reservoir degrees of freedom to be traced over, nor to any Born-Markov 
approximation lO . 



A. Correspondence with classical rate equations 



As master equations constitute a way of introducing classical concepts into a quantum mechanical context, is seems justified 
to hope that the above formalism is able to reproduce the corresponding rate equations for the decay process. Note, however, 
that there is no a priori guarantee that rate equations are relevant here; we are dealing with a quantum mechanical process in 
which not only probabilities, but also coherence effects may play a role. To the extent that the population of the various states 
involved do follow rate equations, the populations should fulfill 

Acs = -rPres (13) 

P P = +r p P ICS , (14) 
which, with initial conditions P ICS (t = 0) = 1, P p (t = 0) = Vp, have the solutions 

P cs = e- rt , P p = ^(l- e - rt ) . (15) 

Here P res (t) is the population of the initial meta-stable state, and P p (t) is the population of the p-th decay product, which is 
assumed to be stable. In order for total population to be conserved, the partial rates r p should sum up to the total rate; 

5> p = r. (16) 

p 

The (partial) rates r( p ) are directly related to the (partial) widths as they simply differ by a factor h; the terms "rates" and 
"widths" are often used interchangeably. 

We may use the complex spectrum of Hqap to identify resonances HMdl]. We will take an A^-particle resonance state ipres^ 
with 

H CAP |e ) }=^ ) ie ) > 
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to be our initial state, p(t = 0) = lores' } (V'rcs'' | ■ We have also assumed that the system is stable after emission of a single particle, 
i.e., the state ipres 1 is such that asymptotically only one particle may be found. As the particle belonging to the single particle 
continuum escapes, it is absorbed, and the (N — l)-particle state pn-i is reconstructed via the source term, cf. Eqs. ( 191 1 Ob . We 

will denote the eigen states of the (N — l)-particle system by ip p N , As the only unbound particle is removed, only bound 
(N — l)-particle states may be populated. (This may be seen from Eqs. (19110b .) Since these states have a negligible overlap with 
the CAP, they all have real eigen energies, and they are all orthogonal to each other; 

ffcAPK- 1 '} = ff|<- l ))= e f- l )K- l )> 

Irnef*- 1 ' =0 Vp (17) 
( ^-i) | ^-x )) = §pg 

The eigen energy of the initial state, however, is complex, and its imaginary part is negative 

Im E W<0, defining = -Im^' . (18) 

The population of the initial state P res is found from the density matrix as (ipr^ \p\ip^el?) = (ipTes\pN\wei?)- With this, 
Eq. (fTTT ) provides 

ihP -£ {N) P -P (- {N) Y = -He 1 P 

nil res — t res ± ros res I t-ros I z,tt res-' res 

P - --e 1 P 

I.e., Eq. (13[ is reproduced with the familiar identification r = 2/K e{ cs , or in terms of a width: 

r = 2 e rcs 

The explicit solution of Eq. ( fTTT i with the proper initial condition is 

p N (t)=e- rt \4£){4» a >\ . (19) 
The population of the p-th (N — l)-particle eigen state Lp\^~ x ^ is 

P p (t) = Tr [H^W^Mt)} - H^W-il^) • 
With this, Eqs. ( 191121 ) and the fact that all the involved eigen energies are real, Eq. ( fT7] i, we have 

l nP p = 4^ i )p p -p p (4 w - i ))* + l ft(<- 1) |5K]i<- l )> # 

P p = (^f- 1 )|5[p Ar ]|<- 1 ))=e- rt (<- 1 )|5[|^ ) )(^ ) l] K" 1} ) = 

kl 

We have here used Eq. ( fT2l . Eq. ( TT4| is reproduced by identifying r p with (ip p N ^ \S[\ipreP)('4'ies > WlPp* )■ Thus, the partial 
width for the p-th decay product is B45I1 : 

r p = 2 J2(x k \r\xi) H^l^) (^UH"^) • (20) 

kl 

Partial widths defined as in Eq. (f20]l will manifestly sum up to the total width, cf. Eq. ( TToT l. This is ensured by the fact that 
the Lindblad equation is trace conserving and that the relevant (N — l)-particle eigen states are orthogonal, cf. Eq. ( fTTI i. Also 
the reality and positivity of the partial widths is ensured by the the fact that they are derived from an equation of Lindblad form. 
This is also clearly seen from a diagonal representation, i.e. a grid representation, of the CAP, Eq. ( ITUl ): 



(21) 



Conceptually, Eq. d20t offers a rather appealing way of calculating partial widths as there exists such a clear correspondence 
between the presented formalism and the classical rate equations for the populations. 
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B. Loss of coherence 

As information of the escaping particle is removed by the CAP, also the description of the system in terms of a pure state is 
lost. This loss of coherence may be quantified by the von Neumann entropy, S(p) = — Tr[plnp], or the purity q(p) = Trp 2 . 
For our analysis we find that the latter quantity, which is unity for a pure state and decreasing with the degree of "mixedness", 
will be the most convenient one. With our density operator expressed as p = px + Pn-i, where p^ is provided by Eq. < TT~9b and 
Pn-i is expressed in terms of the (stable) eigen states, 

PN - 1 =Y J Pir i \tU N) )(^ N - 1) \ , (22) 

r,s 

the purity may be expressed as 

? (p) = e -^ + ^| p (^-l)|2 . (23) 

rs 

The coefficients Prs' - may be obtained from Eq. (O e.g. by means of the Laplace transform. With the proper initial conditions, 
the solution is 

= {e -iAe r .t/n_ e -rt/n ) with (24) 

1 — ii\e rs 

Ae rs = e^-V - ei N -V e E and 

Krs = 2^( Xfc |r| Xi )(^- 1) |c^S ) )(^ ) l4l^- 1) ) ■ 

M 

The diagonal elements, or populations, p rr assumes the simpler form provided in Eq. (15[ (n rr = T r ). Thus, the purity, as a 
function of time, is provided by 

rs rs rs rs 

where equality applies when there is a single final state accessible only. 

III. PARTIAL WIDTHS IN THE CONTEXT OF EXTERIOR COMPLEX SCALING 

As mentioned, there exist several other techniques for obtaining normalizable representations of resonance wave functions. A 
rather straight forward generalization of the above formalism can accommodate for exterior complex scaling. Exterior complex 
scaling is introduced by modifying the spatial variable r such that the radial distance from the origin, r, is rotated by an angle 9 
into the first quadrant of the complex plane for positions beyond a certain radial distance Ro, i.e. r — > R(r) such that 

l#( r )l = J n , iBi D , ' "n ■ (26) 
1 Wl \ R + e l6 {r - Rq), r > R 

Just as in the case of a CAP, stable states, i.e. localized states, are virtually unaffected by the scaling as long as the "unsealed 
region" is large enough, i.e. as long as Rq is larger than the extension of all populated bound states. With an adequate imple- 
mentation of ECS, distinguishing resonance states from other continuum states is rather straight forward. "Ordinary" continuum 
states are typically rotated by the angle 26 into the fourth quadrant of the complex plane from their respective thresholds, whereas 
resonance states are virtually ^-independent and hence easily separated from the other continuum states given large enough 6. 

By extending the formalism for CAPs more or less directly to Hamiltonians modified by ECS, we may arrive at a formula for 
partial widths valid for this case. If the unsealed Hamiltonian is represented in terms of single particle orbitals, Xu ECS will 
introduce an anti-Hermitian term in the Hamiltonian which reads 

H ah = £ hljcld + l - ]T V^ rs clc\c s c r (27) 

kl pqrs 



where 



hi, = / d 3 r Oa(r))* 

Jr>R 



h 2 

sin(260— V 2 + ImVx(R(r)) 
2m 



Xi(r) and (28) 
V^ trs = I d 3 rdV( Xp (r))*( X? (r'))*[Im^(i?(r),i?(r'))]x r (r) Xs (rO . (29) 



r>R Jr'>R 
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Here V\ (r) is the local one-particle potential and V 2 (r, r') is the interaction. Of course, the field operators may still also refer to 
grid-points. In such a context, care must be taken when constructing the coefficients corresponding to the kinetic energy operator 
from some (high-order) finite difference scheme such that the cusp condition at |r| = Rq is fulfilled. 

Reasoning completely analogous to what was performed in the Sec.|II]leads to the following equation for the evolution of the 
n-particle system: 

ihp n = [H h , Pn ] -i{H ah ,p n } +iHS 1 [p n+1 ) +ihS 2 [p n +2} with (30) 
Si[p n +i) = t h l k l ci p n+1 4 and (31) 

kl 

5 2 [p n+2 ] = ^^2v^ qrs c s c r p n+2 clcl . (32) 

pqrs 

If the initial system is represented by an iV-particle pure state we have the following special cases: 

Pn = with iMf^=H^ and 

ihp N -! = F/9JV-1 - PN-1& + ihS 1 [\V N )($ N \] . 

Arguably, the biggest difference between this approach and the CAP-approach is the fact that Eq. ( f30b features a two-particle 
source term, as opposed to Eq. (O which only includes a one-particle source term. The "flow" is still diagonal, however; there are 
no elements in the total density matrix which correspond to a particle number which is not well defined (unless such states were 
populated initially). If we again are to consider a resonance state which is stable after emission of a particle, the two-particle 
source term is not of crucial importance, however. The initial "stationary" state ^res will have at most one particle in the region 
beyond Rq if Ro is chosen large enough. Since V'rcs' is an eigen state of the non-Hermitian Hamiltonian, this situation will 
remain. Thus, as ^[pw] removes two particles in the region beyond Rq, its contribution vanishes in this situation. 

In the same manner as above, the classical rate equations, Eqs. ( 1131141 1. are reproduced by Eq. ( f30b . The partial rate r p = T p /h, 
i.e. the rate at which the p-th eigen state of the (N — l)-particle system is populated, is now calculated as 

r p = 2^^(^- 1 )|c i |^ ) )« ) l4l^ JV) ) . 03) 

k,l 

where h\ ; is provided by Eq. d28| i. Technically, the only difference between the above formula and Eq. ( f20T > lies in how the 
coefficients (xfe|r|xi), resp. h\ t , are obtained. 



IV. DISCUSSION 



Let us first address the advantages with the formalism presented here. The greatest numerical advantage lies, as in 113011 . in the 
fact that rather than analyzing an unbound A -particle wave function, we may now obtain the relevant information by studying a 
localized (N — l)-particle system - once the resonance state is obtained. This is a considerable simplification for two reasons. 
First, the "curse of dimensionality" is reduced. Second, partial widths are obtained without explicit reference to scattering states; 
the finite set of bound states suffices once the initial resonance wave function has been determined. The fact that the population 
dynamics is dictated by an equation of Lindblad form ensures that all partial widths sum up to the total width. 

For the predictions of Eq. ( 120b or Eq. ( l33l to be reliable, the initial resonance state must be well represented. The complex 
resonance energy should remain constant as the absorber region is moved outwards. Moreover, the change induced by a reduction 
in the strength of the CAP or, for ECS, an increase in the scaling angle 6 should be minimal. It must be checked in any 
implementation to what extent these criteria are met. Clearly, the partial widths T p are subject to the same criteria for invariance 
as the total width V. 

Some issues remain unresolved. Possibly the most obvious one is that no formalism has been presented for (uniform) complex 
scaling, which is imposed simply by multiplying the position vector by a complex phase, r — > re l6 . A generalization of the 
present formulae which also accommodates for this form of non-Hermicity is desirable. This is not trivially obtained, however. 
Another rather unsatisfactory feature is revealed if we consider an initial state which is not a single resonant state. Supposing 
our initial state may be written as the super position v&W = citp^ + C2</9 2 '■ If 2 are eig en states, it would seem natural 
to interpret the incoherent sum \c\\ 2 + \c 2 \ 2 as the total population of the initial A-particle system. However, as the states are 
eigen states of a non-Hermitian Hamiltonian, they will in general not be orthogonal, and the sum |ci| 2 + \c 2 \ 2 will not provide 
the trace of p^. Thus, it is not straight forward to interpret such a state in terms of populations. 

Another interesting question arises if we omit the assumption that the system is stable after emission of a single particle, 
i.e. the initial resonance state does not exclusively belong to the single-particle continuum. Can we still use these concepts to 
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describe unstable systems which decay into sub-systems which themselves may be unstable? A straight forward application 
would not provide reasonable results for the same reason as above; the wave functions of the intermediate decay products, 
<Pp N , would in general not be orthogonal (Im e p N ^ < 0). 

All these issues call for a generalization of the formalism presented here. It seems that such a generalization should be based 
on a bi-orthogonal basis representation rather than an orthogonal one 14411 . 

V. CONCLUSION 

Formulae for calculating partial widths for decay-processes in which the remaining system is stable after emission of one 
particle was derived. The formalism combined complex absorbing potentials and exterior complex scaling, respectively, with 
second quantization and the Lindblad equation. The partial widths was derived in a consistent manner which ensured that they 
sum up to the total width of the original meta-stable state. The proposed way of calculating partial widths should be rather 
easily implemented once a good representation of the initial resonance state is obtained, and it is also believed to be numerically 
favorable. 
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